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1 Introduction and summary 



In supersymmetric (susy) theories with sohtons, topological quantum num- 
bers appear as central extensions of the susy algebra. In the topological sec- 
tors there is then typically a lower bound to the mass spectrum determined by 
the central charge, and BPS (Bogomolnyi-Prasad-Sommerfield) states, which 
saturate this lower bound, form shortened multiplets when compared to the 
usual massive multiplets. They are of particular interest because this "mul- 
tiplet shortening" ties the mass of these states to their topological quantum 
numbers such that the BPS saturation is protected from quantum corrections 
[1]. (This is sometimes overstated as implying that there are no quantum 
corrections to the classical mass spectrum at all, while it rather means that 
such quantum corrections have to affect mass and central charge by equal 
amounts.) 

The simplest example of a susy theory with solitons, the H-l dimensional 
susy kink, seemed to be an exception in that the counterparts of short and 
long multiplets have an equal number of states (namely two) [1] according 
to standard representation theory of susy. Nevertheless, most of the explicit 
calculations found neither nontrivial corrections to the mass [2] nor to the 
central charge [3]. 

A couple of years ago, this issue was reopened when two of the present 
authors found [4] that the simple energy-momentum cutoff used explicitly or 
implicitly in most of the calculations that obtained a null result was incon- 
sistent with the integrability of the bosonic sine-Gordon model [5] . Careful 
calculations using mode number cutoff such that boundary energy contribu- 
tions arc avoided subsequently established a nonzero result for the quantum 
corrections to the mass [6, 7, 8, 9] that agreed in fact with an older result 
by Schonfeld [10] and which has also been reproduced by different methods 
[11, 12, 13]. On the other hand, it was established by Shifman et al. [14], 
using a susy-preserving higher-derivative regularization method, that there is 
also an anomalous contribution to the central charge which still leads to BPS 
saturation at the quantum level^, and which was subsequently explained by 
the possibility of "super-short" single-state supermultiplets in 1-1-1 dimen- 
sions [16, 8] which have no definite fermion number. 

In Ref. [17] we have recently shown that these results are most elegantly 

^ A novel derivation of the central charge anomaly using superspace methods to evaluate 
the Jacobian in the path integral with heat-kernel methods can be found in [15]. 



1 



derived by employing dimensional regularization through susy-preserving di- 
mensional reduction from a higher-dimensional model. The correct quantum 
corrections to the mass of the susy kink are obtained [18] without having to 
deal with energy located at boundaries introduced in other methods, and the 
anomalous contribution to the central charge can be obtained from correc- 
tions to the momentum operator in the extra dimensions, which in the case 
of a kink background leave a finite remainder in the limit of 2 dimensions 
[17]. 

In this paper we consider the Abrikosov-Nielsen-Oleson [19, 20, 21, 22] 
vortex solution of the abelian Higgs model in 2-1-1 dimensions which has a 
supersymmetric extension [23, 24] (see also [25, 26]) such that classically the 
Bogomolnyi bound [27] is saturated. We employ our variant of dimensional 
regularization to the N = 2 vortex by dimensionally reducing the = 1 
abelian Higgs model in 3 + 1 dimensions. We confirm the results of [23, 
28, 29] that in a particular gauge (background-covariant Feynman-'t Hooft) 
the sums over zero-point energies of fluctuations in the vortex background 
cancel completely, but contrary to [23, 28] we find a nonvanishing quantum 
correction to the vortex mass coming from a finite renormalization of the 
expectation value of the Higgs field in this gauge [30, 29]. In contrast to 
[23], where a null result for the quantum corrections to the central charge 
was stated, we show that the central charge receives also a net nonvanishing 
quantum correction, namely from a nontrivial phase in the fiuctuations of 
the Higgs field in the vortex background, which contributes to the central 
charge even though the latter is a surface term that can be evaluated far 
away from the vortex. The correction to the central charge exactly matches 
the correction to the mass of the vortex. 

In Ref . [28] , it was claimed that the usual multiplet shortening arguments 
in favor of BPS saturation would not be applicable to the N = 2 vortex since 
in the vortex background there would be two rather than one fermionic zero 
modes [31], leading to two short multiplets which have the same number of 
states as one long multiplet.^ We show however that the extra zero mode 
postulated in [28] has to be discarded because its gaugino component is 
singular, and that only after doing so there is agreement with the results from 
index theorems [32, 31, 33]. For this reason, standard multiplet shortening 
arguments do apply, explaining the BPS saturation at the quantum level that 

^Incidentally, Refs. [28, 31] considered supersymmetric Maxwell-Chern-Simons theory, 
which contains the supersymmetric abelian Higgs model as a special case. 
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we observe in our explicit one-loop calculations. 

2 The vortex in 3 < D < 4 dimensions 

The N = 2 susy vortex in 2+1 dimensions is the solitonic (finite-energy) 
solution of the abelian Higgs model which can be obtained by dimensional 
reduction from a 3+1-dimensional N — 1 model. We shall use the latter 
for the purpose of dimensional regularization of the 2+1-dimensional model 
by susy-preserving dimensional reduction from 3+1 dimensions (where the 
vortex has infinite mass but finite energy-density). 

2.1 The model 

The superspace action for the vortex in terms of 3+1-dimensional superfields 
contains an = 1 abelian vector multiplet and an = 1 scalar multi- 
plet, coupled as usual, together with a Fayet-Iliopoulos term but without 
superpotential, 

jC^ J d'^eww^^ J d'^e^e^^^ + K J d^ev. (i) 

In terms of 2-component spinors in 3+1 dimensions, the action reads^ 

>C = -\Fl. + t^-<0^,x' + \D' + {n-e\ct>\')D 

-\D^ct>\^ + riKpD^^^ + + V2e [0*XaV'" + ^XdV^I , (2) 

where = — ieA^ when acting on and ip, and F^^, = d^Ay — d^A^. 
Elimination of the auxihary field D yields the scalar potential V — \D'^ — 

ie2(|0|2 - t;2)2 with v'^ = k/c. 

In 3+1 dimensions, this model has a chiral anomaly, and in order to 
cancel the chiral U(l) anomaly, additional scalar multiplets would be needed 
such that the sum over charges vanishes, Ylii^i = 0- present paper 

we shall consider only the dimensional reduction of the minimal model (2), 

^Our conventions are r]'^" = (— 1, +1, +1, +1), x" = e'^^Xff and x" = e^'^X/j with 
e"^ = ea/3 = — e"'^ = — e^^ and e^^ = +1. In particular we have •ipa = (V'a)* but 
■ip" = —{ip")*. Furthermore, = {—l,a) with the usual representation for the Pauli 
matrices a. 
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postponing a discussion of anomaly-free 3+1 dimensional vortices [34, 25] to 
a forthcoming work [35]. 

In 2+1 dimensions, dimensional reduction gives an = 2 model involv- 
ing, in the notation of [28], a real scalar N = A3 and two complex (Dirac) 
spinors ip = (^"), x = (x")- 

Completing squares in the bosonic part of the classical Hamiltonian den- 
sity one finds the Bogomolnyi equations and the central charge 

= ^\Dkcl> + ie« + I {F12 + e(|0r - v'))' 

+ -v^€kiFki - idk{eki(j)*Di(j)) (3) 

where k, I are the spatial indices in 2+1 dimensions. The classical central 
charge reads 

j (fx ekidk {ev'^Ai - i(t>*Di4>) , (4) 

where asymptotically Di(p tends to zero exponentially fast. Classically, BPS 
saturation E — \Z\ — 27rv'^n holds when the BPS equations {Di ± iD2)<p = 
D±(f) — and F12 ± e(|0p — v'^) — are satisfied, where the upper and lower 
sign corresponds to vortex and antivortex, respectively. The vortex solution 
with winding number n is given by 

0v = e-7(r), eAY = -ze'^^^^i^, Al = AY±tA^ (5) 

r 

where /'(r) = -/(r) and a'(r) = re'^{f{ry — v^) with boundary conditions 
[22] 

a(r — i> 00) = 0, /(r 00) = v, 
a(r ^ 0) = n + 0{r^), f{r ^ 0) oc r" + 0{r"+'^). (6) 



2.2 Fluctuation equations 

For the calculation of quantum corrections to a vortex solution we decompose 
(p into a classical background part 0v and a quantum part rj. Similarly, A^ 
is decomposed as A^ + a^, where only with fj, — 1,2 is nonvanishing. 
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We use a background [36] gauge fixing term which is quadratic in the 
quantum gauge fields, 

c,.^ = -^i^a' - ^eii<t>YV* - Mf- (7) 

The corresponding Faddeev-Popov Lagrangian reads 

>Cghost -h{dl- {2 + 0vry* + (t>wr]]) c . (8) 

The fiuctuation equations in 2+1 dimensions have been given in [23, 28] 
for the choice ^ = 1 (Feynman-'t Hooft gauge) which leads to important 
simplifications. We shall mostly use this gauge choice when considering fluc- 
tuations in the solitonic background, but will carry out renormalization in 
the trivial vacuum for general ^ to exhibit some of the intermediate gauge 
dependences. 

Because we are going to consider dimensional regularization by dimen- 
sional reduction from the 3+1 dimensional model, we shall need the form of 
the fluctuation equations with derivatives in the direction included. (This 
one trivial extra dimension will eventually be turned into e — > dimensions.) 

In the 't Hooft-Feynman gauge, the part of the bosonic action quadratic 
in the quantum fields reads 

4ol = -^(w-e^i0vix-i^.^r-^'(3i0vr-.^)ir^r 

-2iea^ [r)*Dl4>w - r?(i^>v)1 ■ (9) 

In the trivial vacuum, which corresponds to (j)Y ^ v and — > 0, the last 
term vanishes, but in the solitonic vacuum it couples the linearized field 
equations for the fiuctuations B = {r],a^/y/2) with a+ = ai + ia2 to each 
other according to (A; = 1, 2) 

^"^^-"^tJ^-y -^^/2e(D_0v)* -9^ + 2e^|</.v| V ^ ^ 

The quartet (03, ao, b, c) with b, c the Faddeev-Popov ghost fields has diagonal 
field equations at the linearized level 

(dl - 2e2|0v|')g = 0, Q= (as, ao, b, c). (11) 
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For the fermionic fluctuations, which we group as U — (^^i V 
the hnearized fleld equations read 

LU = i{dt + d3)V, LV = i{dt - ds)U, (12) 

with 



-V^e(t)*y id_ y ' V A/2e0v id+ )' ^^^^ 

Iteration shows that U satisfies the same second order equations as the 
bosonic fluctuations 

L^LU={dl-df)U, L^LB^{dl-dl)B (14) 
LL'^V = {dl - d',)V, (15) 

with L^L given by (10), whereas V is governed by a diagonal equation with 

(In deriving these fluctuation equations we used the BPS equations through- 
out.) 



2.3 Renormalization 

At the classical level, the energy and central charge of vortices are multiples 
of 27rf ^ with v"^ = n/e. Renormalization of tadpoles, even when only by finite 
amounts, will therefore contribute directly to the quantum mass and central 
charge of the N — 2 vortex, a fact that has been overlooked in the original 
hterature [23, 28] on quantum corrections to the N — 2 vortex.^ 

Adopting a "minimal" renormalization scheme where the scalar wave 
function renormalization constant Z^p = 1, the renormalization of v"^ is fixed 
by the requirement of vanishing tadpoles in the trivial sector of the 2+1 di- 
mensional model. The calculation can be conveniently performed by using 
dimensional regularization of the 3-1-1 dimensional N = 1 model. For the cal- 
culation of the tadpoles we decompose = v+r] = v+{a+ip) / a/2, where a is 
the Higgs field and p the would-be Goldstone boson. The gauge fixing term 

"^The nontrivial renormalization of k/c has however been included in [30, 29]. 
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(7) avoids mixed a^-p propagators, but there are mixed X'i^ propagators, 
which can be diagonahzed by introducing new spinors s = (?/' + ix)/\/2 and 
d = {ijj — ix)/V^ with mass terms m(sQ,s" ^ da<P) + /i.e., where m = \/2ev. 

The part of the interaction Lagrangian which is relevant for a tadpoles 
to one-loop order is given by 



em 



C-tadpoies = <x.r + Xc.r)(y-^{^' + p') a-em{a'^Hbc-5v') a, (17) 



where b and c are the Faddeev-Popov fields. 

The one-loop contributions to the a tadpole thus read 



s,d 




'^0.1,2.3 





b,c 



5v^ 



-em) X 



(18) 



3 1 1 1 1 

{-2trl2/(m) + -/(m) + -/(C^m) + [3/(m) + ^I{^^m)] - ^I{^^m) - 6v^}, 



where 
7(m) 



mi+^ r(- 



(27r)3+^ A;2 + m^ ~ (47r)i+^/2 Y{-\ 



^--^ + 0{e). (19) 



Requiring that the sum of tadpole diagrams (18) vanishes fixes Sv'^, 



Sv'^ = ^ (^I{m) + I{C^m)^ 



D=3 



-m. 



(20) 



Because in dimensional regularization there are no poles in odd dimensions at 
the one-loop level, the result for Sv^ is finite, but it is nonvanishing. Because 
the classical mass of the vortex is My — 27rv'^ — 7rm^/e^, the counterterm 
Sv"^ is the only one that is of importance to the one-loop corrections to My. 
Since 5f ^ is gauge-parameter dependent, the remaining contributions to mass 
and central charge must be gauge dependent, too, so that the final result is 
gauge independent. 
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3 Quantum corrections to mass and central 
charge 

3.1 Mass 

At the one-loop level, the quantum mass of a solitonic state is given by 

M = Mci + ^ ^ cu^os -\Y1 ^^^'^ + (21) 

where is the classical mass expressed in terms of renormalized parame- 
ters, 5M represents the effects of the counter-terms to these renormahzed pa- 
rameters, and the sums are over zero-point energies in the soliton background 
(the zero-point energies in the trivial vacuum, which one should subtract in 
principle, cancel in a susy theory). 

In the ^ = 1 gauge the sum over zero-point energies is formally 

2 ^ ~ 2 ^ '^ferm = ^ + ^ — ^ — ^ 

= (22) 

where the quartet {as,ao,b,c) cancels separately. (Note that in (22) all fre- 
quencies appear twice because all fields are complex.) Using dimensional 
regularization as developed in [18], these sums can be made well defined by 
replacing all eigen frequencies LOk in 2-1-1 dimensions by Uk^e — {^l + i'^Y^'^ 
where £ are the extra momenta, and integrating over £. The spectral den- 
sities are nontrivial only with respect to the 2-dimensional momenta k and 
the former arc not modified by dimensional regularization. 

In [23, 28] it has been shown that L'^L and LL"!", which govern U and V, 
respectively, are isospectral up to zero modes. In dimensional regularization, 
where the zero-mode contributions continue to give zero because scaleless 
integrals vanish, one can therefore conclude that in the ^ = 1 gauge there 
is a complete cancellation of the sums over zero-point energies. All that 
remains is the finite renormalization of in that gauge: 

E = 2n\n\iv' + 5v\^,) = 2n\n\iv' - £) ^ \n\ - ^) (23) 

(In gauges other than ^ = 1 the fiuctuation equations for the B fields, i.e. 
T), a+, no longer match those of the U fermions.) 
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This result agrees with [29] , where however a careful analysis of boundary 
conditions in the heat-kernel approach was needed because the vortex had 
to be put in a box to discretize the spectrum. In dimensional regularization 
one does not need to put the system in a box, and as a consequence there is 
no need to study the contributions from these artificial boundaries. 

3.2 Central charge 

By starting from the susy algebra in 3+1 dimensions one can derive the 
central charge in 2+1 dimensions as the component T03 of 



where is the susy Noether current. 

The antisymmetric part of T^" gives the standard expression for the cen- 
tral charge density, while the symmetric part is a genuine momentum in the 
extra dimension: 



(A similar decomposition is valid for the kink [17].) 

Z corresponds to the classical expression for the central charge. Being a 
surface term, its quantum corrections can be evaluated at infinity: 



with (i — ev^Ai — i(f)Wi(f) and Vq — v"^ + 5v^. 

Expanding in quantum fields (j) = (py + r], A — + a and using that the 
classical fields (py ve^"'^, A^ — > n/e, Djcpy ^ as r — > 00, we obtain to 
one-loop order 



(24) 




(25) 




(26) 



(Z) = 2Tinvl -i de + 77t)(DV _ ,eae){(py + r,)) 




Jo 



= 27rn{vQ - {r]^r])\r^oo} - i dO {{ij^'der]) - ie(j)y {aefj) - ie(f)y {aer]^)} 




where we have used {ri{r^ 00)) — > (which determines 5v'^), J^^d6{ao) — 0, 
and (rj'^rjao) — 0{h'^). 
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The first contribution, Z^-, can be easily evaluated for arbitrary gauge 
parameter ^, yielding 

Za = 27rn{t;^ - i((aa) + (pp))|^^oo} 

= 27rnK-^[/(m) + /(e^m)]} 

= 1'Kn{vl - bv^) = l-Kuv^. (28) 

If this was all, the quantum corrections to Z would cancel, just as in the 
naive calculation of Z in the susy kink [3, 4]. 

The second contribution in (27), however, docs not vanish when taking 
the limit r oo. This contribution is simplest in the ^ = 1 gauge, where 
the T] and fluctuations are governed by the fluctuation equations (10). In 
the limit r ^ oo one has |0v| ~^ v and £)_0v exponentially. This 
eliminates the contributions from {aorj). However, D^, which governs the 
77 fluctuations, contains long-range contributions from the vector potential. 
Making a separation of variables in r and 9 one flnds that asymptotically 
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\D^,v\'^\drV\' + -,\(do-in)rj\' (29) 

so that the r] fluctuations have an extra phase factor e*"^ compared to those 
in the trivial vacuum. We thus have, in the ^ = 1 gauge, 

Zh = -i 



f-27r 



^0 

i / d9 {r]'^der]) = 27rn (ry^ry) = 27rn 5v^ . (30) 

JO 



This is exactly the result for the one-loop correction to the mass of the vortex 
in eq. (23), implying saturation of the BPS bound provided that there are 
now no anomalous contributions to the central charge operator as there are 
in the case in the = 1 susy kink [17]. 

In dimensional regularization by dimensional reduction from a higher- 
dimensional model such anomalous contributions to the central charge op- 
erator come from a flnite remainder of the extra momentum operator [17]. 
The latter is given by [28] 

Zc = (Ps) = j d'x {FoiFs, + {Do^yDs4> + {Ds<f>yDo<f> 

-iX^odsX - ii^^oDaip) . (31) 
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Inserting mode expansions for the quantum fields one immediately finds that 

the bosonic contributions vanish because of symmetry in the extra trivial 
dimension. However, this is not the case for the fermionic fields, which have 
a mode expansion of the form 

where we have not written out explicitly the zero- modes (for which uP" — (?\ 
The fermionic contribution to reads, schematically, 

where uj = y/cok + P'-, so that the i integral does give a nonvanishing result. 
However, the ^-integration over the mode functions u\^2^]'^^ and fi,2(fc;a;) 
produces their spectral densities, which cancel up to zero-mode contributions 
as we have seen above^, and zero- mode contributions only produce scaleless 
integrals which vanish in dimensional regularization. Hence, Z^. — ^ and 
\Z\ = \Za + Zb\ — E, so that the BPS bound is saturated at the (one- loop) 
quantum level. 

3.3 Fermionic zero modes and multiplet shortening 

Massive representations of the Poincare supersymmetry algebra for which 
the absolute value of the central charge equals the energy, i.e. when the 
BPS bound is saturated, contain as many states as massless representations, 
which is half of that of massive representations for which the BPS bound is 
not saturated. These results also apply in 2-1-1 dimensions for the N — 2 
super- Poincare algebra [28]. 

A particular multiplet of states is obtained by taking the vortex solution, 
and acting on it with the susy generators of the N = 2 susy algebra, which 
contains two complex charges Q+, Q-, and their hermitian conjugates {Q+Y 
and {Q-Y- One of these charges, annihilates the vortex solution, while 

^An explicit calculation which confirms these cancellations will be presented in [35]. 
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the other one, Q+, is to hnear order in quantum operators proportional to 
the annihilation operator of a fermionic zero mode. 

However, if there indeed is a second fermionic zero mode in the model as 
postulated in [28]^, in second quantization it would be present in the mode 
expansion of the fermionic quartet U and V, 




non-zero modes. (34) 



As a result, there would then be a quartet of BPS states 

\v), a\\v), a\\v), a\a\\v) (35) 

comprising two short multiplets of = 2 susy, which are degenerate and 
together have as many states as one long multiplet without BPS saturation. 
As stressed in [28], the standard argument for stability of BPS saturation 
under quantum corrections from multiplet shortening [1] thus would not be 
applicable. 

However, we shall now show that there is in fact only a single fermionic 
zero mode in a vortex background with winding number n = 1. To this end, 
we first observe that the zero modes must lie in U, because V is governed by 
the operator LL^ of Eq. (16), whose only zero mode solution is Vq = 0. A 

zero mode for U must satisfy LU = 0, and to analyse this equation we follow 
[28] and set ^^{x,y) = —ie'^^^~2~^"'^^u{r) and — e^^^~^2'>^d{r). The equation 
LU = reduces then to 



dr 



V2ef dr + '-^ 




0, (36) 



where / = f(r) and a — a{r) satisfy /' = ^/ and a' — re^{p — e^). Iterating 
this equation yields 



1. 



+ -dr - ^^^^ - 2e'p )- = Q. (37) 



u 



^In the literature one can in fact find two different conventions for indicating the number 
of fermionic zero modes. Like Refs. [28, 33] we only count the number of zero modes in 
the fermionic quartet {U, V) and not additionally those in the corresponding conjugated 
fields (C/t,yl'). One zero mode in this way of counting then corresponds to a pair of 
creation/annihilation operators. Alternatively one may count the zero modes in both 
{U,V) and {U\V'^) and thus ascribe one zero mode to each creation or annihilation 
operator. The latter way of counting is employed for instance in Ref. [26]. 
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Given a solution for u, the corresponding solution for d follows from LU — 0. 

For given j, this equation has two independent solutions, a linear combi- 
nation of which yields solutions which decrease exponentially fast as r — » cxd. 
Hence, both solutions should be regular at r = 0. For j 7^ |, one has, using 
/(r ^ 0) ~ r", 

jj;^ r^u^ r^'idr^-^ + C2r-^^-^^) for r ^ (38) 

which selects for n = 1 only j — — |. This solution is the zero mode that 
is obtained by acting with on the background solutions, which gives 
^lJ^ = -iD„0v/V2 = -iV2f, = F12 = -e{f - v^). For j = i, one finds 
for n = 1 near r = 

^jj^ ^Ci(x + iy) + C2(x + iy)\nr . (39) 

For large r, ijji ~ e~""'e*^, as follows from (37). This solution corresponds to 
the second fermionic zero mode postulated in Ref. [28]. 

However, while (39) is regular at the origin, the associated gaugino com- 
ponent is not: (36) implies that 

X' - C2—, (40) 
r 

so this solution has to be discarded when C2 7^ 0. 

Similarly, one can show that for winding number n > 1 regularity of 
the gaugino component generically requires that J < — | so that the correct 
quantization condition for normalizable fermionic zero modes is —n + ^ < 
j < — |. Hence, there are n independent fermionic zero modes, not 2n as 
concluded in [28]. It is in fact only the former value that agrees with the 
results [31, 33] obtained from the index theorem [32]. 

As has been proved rigorously in [37], in the bosonic sector there are 2n 
zero modes, which are related to the above n independent fermionic zero 
modes by supersymmetry. In the R^=i background gauge c?^a^ — ie{4>yr]* — 
4>yV) = 0, the bosonic zero modes satisfy a set of equations completely equiv- 
alent to those for the fermionic zero modes [31]. But the linearly dependent 
solutions (q) and ^(q) correspond to linearly independent solutions for the 
bosonic zero modes a and 77.^ In particular, for n = 1, the j = — | solu- 
tion {ip^ = —iu{r),x^ — d{r)) with real u{r) and d{r) corresponds to the 

^For an analogous case see eq. (3.8) of Ref. [38]. 
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bosonic zero mode r]{r) — —iu{r), (01,02) = {■\/2d{r) , 0) , while multiplying 
the fermionic solution by i corresponds to the bosonic zero mode ri{r) = u{r), 
(oi, 02) = (0, V2d(r)), which is evidently linearly independent of the former. 
For both solutions the R^=i gauge condition is satisfied due to the lower com- 
ponent of the field equation (36). Conversely, one can start from the classical 
vortex solution and find two independent bosonic zero modes by consider- 
ing their derivatives with respect to the x and y coordinates. Performing a 
gauge transformation to satisfy the R^=i gauge condition leads one back to 
the above solutions. This additionally confirms that the above analysis has 
identified all fermionic zero modes in the quartet {U,V). 

We thus conclude that for the basic vortex (winding number n — 1) there 
is exactly one fermionic zero mode (corresponding to one pair of fermionic 
creation/annihilation operators) and this gives rise to a single short multiplet 
at the quantum level. Standard multiplet shortening arguments therefore do 
apply and explain the preservation of BPS saturation that we verified at 
one-loop order. 
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